Field-tunable stochasticity in the magnetization reversal of a cylindrical nanomagnet 
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The nature of magnetization reversal in an isolated cylindrical nanomagnet has been studied 
employing time-resolved magnetoresistance measurement. We find that the reversal mode is highly 
stochastic, occurring either by multimode or single-step switching. Intriguingly, the stochasticity 
was found to depend on the alignment of the driving magnetic field to the long axis of the nanowires, 
where predominantly multimode switching gives way to single-step switching behavior as the field 
direction is rotated from parallel to transverse with respect to the nanowire axis. 



Traditionally, cylindrical nanomagnets have been of 
great interest for high density magnetic storage [ij, but 
very recently the dynamics of domain walls (DW) in 
these systems is also rapidly gaining in importance. A 
driving factor to this is the absence of "Walker break- 
down" , with the DWs acting as "massless" particles hav- 
ing zero kinetic energy However, a generic problem, 
in both planar or cylindrical nanowires, is the stochas- 
ticity associated with the magnetization reversal pro- 
cess. This is manifested in two classes of phenomena: 
first, the stochasticity associated with diffusive and non- 
deterministic motion of the DWs in presence of arti- 
ficial or intrinsic disorder (addressed elsewhere by the 
same authors Q); second, the stochasticity related to 
the nucleation and subsequent propagation of DWs. The 
latter has been investigated in several planar magnetic 
films and lithographically patterned nanowires through, 
for example, size and shape distributions of Barkhausen 
avalanches, or extraordinary Hall effect etc [J]. Asym- 
metry in the mode of magnetization switching process on 
reversal of magnetic field polarity has been investigated 
using magnetic force microscopy (5| , which might as well 
be related to the stochasticity issue (MFM probes only 
a fraction of the whole sample). The stochasticity in 
magnetization reversal in case of cylindrical nanowires, 
however, is relatively unexplored. 

Cylindrical magnets have been the classical template 
for theoretical study of the mode of magnetization re- 
versal and more importantly the angular variation of 
the nucleation field, albeit in the limit of an infinite 
and isotropic cylinder Q. Many experiments on mag- 
netization reversal in nanoscale magnetic systems have 
looked at the angular dependence assuming an infinite 
and isotropic cylinder (or strips) , although finite size and 
anisotropy related effects can be crucial in those cases 
15[- Wegrowe et al. could fit their data on angular de- 
pendence with the curling mode prediction for an infinite 
cylinder assuming an activation volume with aspect ratio 
2 : 1 (Ref. [3l). Moreover, surface anisotropy or struc- 
tural defects can also play crucial role and hence cannot 
be treated within the framework of isotropic magneti- 
zation 16|. The principal objective of this article is to 
explore the effect of angular variation on the stochastic- 
ity of magnetization reversal, rather than the variation 
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FIG. 1: (A) Magnetic field dependence of resistance at dif- 
ferent angles between the applied magnetic field and the 
nanowire axis suggesting that the angle dependence of Hsw 
has 180° periodicity. (B) SEM micrograph of a device. (C) 
The switching field histograms for 6 = Q and 6 — 80°. (D) 
The time dependence of resistance at a static low magnetic 
field (1.5 Oe) for 6 = is shown. 



in switching field (i?s^), where we show that finite size, 
disorder, local anisotropy etc. play crucial roles. 

Experimentally, we have taken a different approach. 
Instead of magnetic hysteresis measurement such as us- 
ing microfabricated SQUID 7] which is limited to low 
temperature, or magneto-optic Kerr effect (MOKE) Q 
which requires significant averaging over many field cy- 
cles to improve the signal to noise ratio and hence unsuit- 
able for probing stochastic switching, we have measured 
time-resolved electrical resistance around Hg^ to track 
the nucleation of individual DWs at different angles be- 
tween the applied magnetic field (H) and the long axis of 
the cylindrical nanomagnet. Electrical transport in mag- 
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netic nano-cylinders have been reported before 1^ 17 1, 
but primarily to explore the switching mechanism, rather 
than its stochasticity. Our approach and results may also 
be relevant to device applications such as DW logic sys- 
tems 18 1 or spiral turn sensors [l^, which require the 
use of orthogonal magnetic fields. 

Nickel nanowires were electrochemically grown inside 
anodic alumina templates. The average diameter of the 
nanowire is « 200 nm, where strong shape anisotropy 
(aspect ratio ^ 200) aligns the magnetization along the 
long axis of the nanowire. Details of the growth pro- 
cess and structural characterization can be found else- 
where (2ol |. Following growth, nanowires were dispersed 
on flat silicon oxide substrates, and electron-beam lithog- 
raphy was used to form Ti/Au contact pads for electrical 
measurements on single nanowire as shown in Fig. [T]B. 
The resistance was measured at room temperature in the 
four-probe geometry using standard ac lock-in technique. 

The magnetoresistance measured at different angles 
between the applied magnetic field and the current or 
the long-axis of the cylinder show that the switching 
field has a minimum at = and a maximum at 
d - ±90° (Fig. [IK). This 180° periodicity indicates 
curling-type magnetization reversal mechanism expected 
for nanowires of diameter 200 nm at least when H is 
parallel to the long axis ^g]. The Hg^ measurement was 
repeated 50 times. At 6 = 0, Hg^ has a distribution of 
12% around the most probable value of 240 Oe, while at 
= 80°, the distribution is slightly narrower, about 8% 
with a most probable switching field of 720 Oe (Fig.[Tp). 
Although such a narrowing for the transverse direction 
has been observed before [7j for Nickel nanowires of diam- 
eter less than 100 nm below 6K, the overall width of Hgw 
distribution is significantly broader in our case, proba- 
bly due to the larger diameter of the nanowire (which 
could introduce inhomogeneties thereby lowering the lo- 
cal anisotropy barrier) and enhanced role of thermal ac- 
tivation at room temperature. 

For magnetoresistance to be a viable probe to mag- 
netization reversal, it is essential to identify and sep- 
arate the contributions from anisotropic magnetoresis- 
tance (AMR) arising due to spin-orbit scattering, and 
that from the DWs. To achieve this we have studied the 
time dependence of nanowire resistance (i?) at H Hgw ■ 
The distance between the voltage probes placed midway 
along the length of the nanowire is ^ 5/im. The nanowire 
is magnetized in one direction along the long axis and 
then a small magnetic field is applied in the opposite 
direction. The time dependence of R at = 1.5 Oe ap- 
plied parallel to the long-axis of the nanomagnet shows 
stochastic switching back and forth between discrete 
multi-level states in R (Fig: [Tp). The jumps describe 
increase in R from its base value Rq (~ A.78ft) by AR 
or 2AR, where AR « 3 mft and vice versa. Since H 
is constant, AMR or Lorentz contributions to R do not 
change, and hence the observed switching between dif- 
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FIG. 2: AT, II) Time dependence of resistance for two mea- 
surement cycles with sweeping magnetic field at 6 — with 
R{0) ~ 4.78 mfi. The magnetization reversal involves succes- 
sive stochastic irreversible jumps in resistance which are in- 
teger multiples of a Rdw Between two jumps, the resistance 
rises weakly due to AMR. BI,II) Same for — 65°. Resistive 
jump of size 2Rdm is rare at high values of 9. CI, II) Time 
dependence of resistance showing irreversible single-step and 
multimode switching at ^ = 80° for two measurement cycles. 
Absence of AMR leads to flat resistance plateaus. 



ferent R states can be attributed to the DWs traveling 
in and out of the region between the voltage probes. In- 
creasing R by AR and 2AR corresponds to the existence 
of one and two DWs, respectively, between the voltage 
probes. The positive correction AR can be quantita- 
tively understood from the Levy-Zhang model of spin- 
mixing inside the DWs [21| and has been discussed in 
detail elsewhere 

To study the dynamics of magnetization reversal, H 
was swept at a rate of 2 Oe/S across Hgw while R was 
measured as function of time with a resolution of 100 
mS. Below Hgr,,, as H is increased from zero after be- 
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ing saturated in a high H of opposite polarity, the re- 
sistance decreases monotonically. Above TJsu, , successive 
irreversible jumps in the resistance are observed at 6 = 
(Fig:[2K). Single-step switching is rare at 9 = 0, and the 
striking feature is that the most probable jump size is 
an integer multiple of 3 mJ7 (Fig: . The time interval 
between two successive jumps in a single cycle is non- 
deterministic and varies from one field cycle to another 
(Fig. [2]4_, B). The resistance during that time interval is 
not constant due to AMR and rises weakly with time (or 
H). When H is applied along the transverse direction 
single step switching events become more frequent than 
multimode switching (Fig: [2]CI) and the AMR not be- 
ing significant in transverse direction, R stays constant 
between successive jumps (Fig: [2]3II). 

The resistance jumps at integral multiple of AR allows 
us to understand the magnetization reversal mechanism 
at ^? = 0. Since AR corresponds to the resistance of a sin- 
gle DW (vortex type since the 9 dependence of Hsw sug- 
gests curling mode reversal) , it can be safely argued that 
the stochastic multiple jumps in resistance correspond to 
multiple nucleation events, where jumps of size 6 mTl in- 
dicate the nucleation of a pair of DWs. Motion of DWs 
far below the nucleation field or magnetization reversal 
due to nucleation of a DW pair is consistent with the 
time dependent solution of Landau-Lifshitz equation 22 1 
and particularly for transition metal nanowires [1^ . The 
plateaus correspond to DW pinning / small-scale dis- 
placement or even rotation of spins (given the weak up- 
ward trend in the resistance plateaus at = due to 
AMR). Presence of multiple vortex DW suggests that 
magnetization reversal at 9 = occurs through localized 
curling mode with multiple nucleation events. Due to the 
simplicity of device geometry, it is an open question as 
to where the DW nucleation starts (possibly at the wire 
ends or at some local inhomogeneties within the wire) 
although they are eventually detected within the volt- 
age probes. Similar jumps in the longitudinal resistance 
have been observed in FePd nanostructure 2^. Multi- 
ple vortex walls have also been observed in planar NiFe 
nanowire of larger width 15| . 

The stochastic distribution of switching time between 
each irreversible jump was measured over hundred such 
switching events. The switching probability was com- 
puted from the integral over the switching time his- 
tograms for different values of 9 (where multimode 
switching is predominant), which were fitted with a 
stretched exponential function P{t) — exp— (i/r)*^, with 
/? varying between 1.3 — 2 (Fig: [3j3). The experimental 
data for ^ = follows the overall feature of the theoriti- 
cal fit although some deviation is observed for low time 
scale possibly due to small number of switching events in 
that region. Similar results (/3 > 1) have been reported 
elsewhere, which use more elaborate experimental tech- 
niques such as micro-SQUID, magnetic force microscopy 
or MOKE |25l428| . The probability of not switching as 
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FIG. 3: (A) The jump size distribution during the irreversible 
switching for (A) 9 = (B) The enhanced-exponential proba- 
bility of not switching with ^ > 1 (/3 = 1.86, 1.34 for 9 = 0, 65° 
respectively) indicates a relaxing multi-valleyed free energy 
landscape. (C) Jump size distribution for 6 = 80° (D) Proba- 
bility of single-step switching as a function of 9. Inset: Switch- 
ing field distribution for single step and multimode reversal 
for 6 = 80° is shown separately. 



shown in Fig. [3j3 characterized by /3 > 1 suggests re- 
versal to involve correlated thermally activated processes 
over a distribution of barrier heights in the free energy 
landscape, rather than the Neel-Brown picture of thermal 
activation over a single energy barrier. 

The distribution of jump size taken over 50 measure- 
ment cycles at different values of 9 shows that for low an- 
gles such as 9 — 0, the most probable jumps are peaked 
around 3 and 6 mil while for higher angle such as at 
9 — 80°, an additional strong peak appears around 8 — 9 
mil with the distribution around 3 mil being significantly 
broadened (the resistive jumps of size 6 mil are rare). It 
is readily observed that in contrast to the preponder- 
ance of multimode switching at lower angle, single-step 
switching is more frequent at higher angles (Fig: |3lC,D). 

The origin of enhanced single-step switching at higher 
values of 9 remains unclear. Recent micromagnetic sim- 
ulations have shown that application of transverse H re- 
sults in the expulsion of the vortex wall leading to mag- 
netization reversal via coherent rotation P]. However, 
the experimental evidence is merely based on the devia- 
tion of Hsw value from that predicted for curling mode 
at high aiigle. No such deviation has been observed for 
NiFe 13, 13] or Ni nanowires at low temperatures 0. 
Wegrowe et al. attributed the deviation from curling 



mode for cylindrical Nickel nanowire (average diameter 
60 nrn) at room temperature to pinning by surface de- 
fects [16[. At higher 9, the value of Hsw predicted by 
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coherent rotation model is smaller than that from curl- 
ing model In our case, the curling mode predic- 
tion for a prolate spheroid [29| (assuming an exchange 
length df^j. — 40 nm and bulk saturation magnetization 
Ms — 0.485 T; the demagnetization factors along the 
major and minor axes are = 0.00015 and = 0.5, 
respectively) gives Hsw = 140 Oe at 6 = 0. The under- 
estimation of Hsw value at = has been reported for 
nanowires with d/d^x > 3 [30| . The curling model over- 
estimates Hsw at 6* = 80° by a factor of two, while, on 
the other hand, Hsw value predicted by coherent rotation 
model 3l| (where Hsw peaks at 6* = 0, ±90°) is nearly 
an order of magnitude smaller compared to that observed 
experimentally at 9 = 80°. Moreover, we have observed 
that single-step switching on the average involves higher 
Hsw compared to multimode switching (Inset, Fig:[5p). 
Therefore single step switching cannot be associated with 
coherent rotation mode. One possible mechanism for sin- 
gle jump reversal is nucleation with very small activation 
volume and the vortex wall sweeping across the entire 
length of the nanowire [s^ . The question is whether 
such a scenario is plausible in presence of transverse mag- 
netic field. Transverse H can modify pinning or depin- 
ning processes of DWs js^]. Recently, huge enhance- 
ment of DW velocity under strong transverse H has been 
reported Nonetheless, our experiments prove the 
stochastic nature of the mode of reversal at high values 
of 9 with both single-step and multiple-step reversal be- 
ing possible. 

The observation that application of transverse H on a 
cylindrical nanomagnet has a higher probability to avoid 
multimode switching could prove to be extremely impor- 
tant for device applications. Large area nucleation pads 
are used for DW injection into a nanowire where multi- 
mode switching due to transverse magnetic field has been 
a major obstacle Bryan et al. observed multimode 
switching in wide planar Permalloy nanowire, due to ap- 
plication of transverse H Q. Instead of planar wires 
of rectangular cross-section, nucleation pads with square 
cross-section (to circumvent lithographic obstacles in fab- 
ricating cylindrical nanowire) can be used to reproduce 
the properties of a cylindrical nanomagnet. Use of softer 
magnets such as permalloy can reduce the switching field 
value in the transverse direction. 

To conclude, we have employed a simple time depen- 
dent magnetoresistance measurement technique to probe 
the stochasticity of nucleation-mediated magnetization 
reversal in a cylindrical nanomagnet. When the applied 
magnetic field is parallel to the long axis of the nano- 
magnet, the magnetization reversal essentially follows the 
localized curling mode with more than one nucleation 
events, whereas when the magnetic field is applied per- 
pendicular to the long axis of the nanomagnet, the mode 
of reversal is strikingly non-deterministic and can either 
follow a single-step switching or a multimode switching 
process. 
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